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ABSTRACT 

The evolution of an accretion disk can be influenced significantly by the deposition 
of mass and angular momentum by an infalling Bondi-Hoyle wind. Such a mass influx 
impacts the long-term behavior of the disk by providing additional sources of viscosity 
and heating. In this paper, we derive and solve the disk equations when these effects are 
taken into account. We present a survey of models with various wind configurations and 
demonstrate that the disk spectrum may then differ substantially from that of a standard 
a-disk. In particular, it is likely that a wind-fed disk has a significant infrared bump 
due to the deposition of energy in its outer region. We apply some of the results of our 
calculations to the Galactic Center black hole candidate Sgr A* and show that if a fossil 
disk is present in this source, it must have a very low viscosity parameter (a <^ 10 -4 ) and 
the Bondi-Hoyle wind must be accreting with a very high specific angular momentum 
to prevent it from circularizing in the inner disk region where its impact would be most 
noticeable. 

Subject headings: accretion, accretion disks, Bondi-Hoyle — black hole physics — Galaxy: 
center — Galaxy: Sgr A* 
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1. Introduction 

1.1. Disk vs. Spherical Accretion 

Compact objects, such as white dwarfs, neutron 
stars, and black holes may accrete either from the am- 
bient medium (generally the Interstellar Medium — 
ISM) or from a companion in a binary configuration. 
Whether or not this accretion is mediated via a disk 
or it constitutes a more or less quasi-spherical infall 
depends critically on the specific angular momentum 
carried by the captured gas. It therefore is not al- 
ways clear whether disk or quasi-spherical accretion 
dominates or if both are present. 

For example, in Roche Lobe overflow, the gas cir- 
cularizes quickly and forms a disk that can extend as 
far out as the binary orbital midpoint (e.g., Lamb & 
Melia 1987; Melia & Lamb 1987). Its temperature is 
generally too low to provide any significant support in 
the vertical direction, and so the disk is relatively flat 
and optically thick, at least until the gas approaches 
the inner regions close to the accreting object. Some 
binary systems are too wide for the companion to fill 
its Roche Lobe, and for these the accretion proceeds 
primarily via the capture of a portion of its wind by 
the compact object. In many circumstances, the lat- 
ter's magnetic field is quite strong and the plasma is 
entrained onto its polar cap regions (e.g., Burnard, 
Klein & Arons 1988). Otherwise, the gas can either 
form a small (highly variable) disk whose size is spec- 
ified by the circularization radius dictated by the spe- 
cific angular momentum, or it may fall directly onto 
the accretor (Taam & Fryxell 1988). 

In active galactic nuclei (AGN), the central, mas- 
sive black hole generally accretes from the ambient 
medium. This can either be the ISM gas surround- 
ing the nucleus or it may be the winds produced by 
stars in the central cluster. Although it is likely that 
Bondi-Hoyle accretion dominates the gaseous motions 
at large radii (out at the so-called accretion radius of 
order 0.1 — I pc), there is ample evidence that much 
of the activity in AGNs is associated with an accre- 
tion disk at smaller radii. It appears, therefore, that 
the quasi-spherical accretion at large radii gradually 
settles into a planar, more or less Kcplcrian structure 
toward the central accretor, i.e., the accretion disk in 
these sources is 'fed' by the quasi-spherical infall. 



1.2. Scope of the Paper 

In this paper, we begin our study of the effects 
on the accretion disk evolution due to the continued 
infall from quasi-spherical accretion. We rederive the 
disk equations taking into account the contribution 
from the wind to the overall disk dissipation, the disk 
temperature, its column density, the accretion rate 
through the disk and the overall disk spectrum. All of 
these are expected to differ substantially between the 
isolated systems fed within a plane from their outer 
radii and those in which a disk is mainly fed from 
above. 

The infalling gas can differ from the disk plasma 
in several respects: (i) its specific angular momentum 
at a given radius may be close to Keplerian (like that 
of the disk), or it may be much less, perhaps even 
zero, (ii) its mass flux may be much larger than that 
through the disk, and (iii) its ram pressure may dif- 
fer significantly from that of the disk material, and 
thereby produce its own dynamical influence on the 
gas closer to the compact object. 

We wish to study the impact of such a mass influx 
on the long-term behavior of the disk. In future work, 
we will couple this analysis with 3D hydrodynamic 
(and possibly magnetohydrodynamic) simulations of 
the ambient infall. This will provide much more re- 
alistic profiles of accreted mass and specific angular 
momentum as functions of radius in the disk. The 
process we discuss may be relevant to galactic nuclei 
as well as stellar mass black holes, neutron stars and 
star formation. Here we will concentrate on the for- 
mer. 

1.3. The Galactic Center 

We shall specifically consider the example of the 
Galactic Center, where the presence of a supermassive 
black hole (coincident with the unique radio source 
Sgr A*) and strong, nearby stellar winds should 
lead to substantial Bondi-Hoyle accretion (e.g., Melia 
1992; Ozernoy 1993; Ruffert & Melia 1994). In addi- 
tion, stars may get tidally disrupted as they approach 
the strong central gravitational field. Depending on 
the strength of the disruption, as much as half of the 
stellar mass may be left as a remnant surrounding 
the black hole (Khohklov & Melia 1995 and references 
cited therein, particularly Rees 1990 and Shlossman, 
Begelman & Frank 1990). This remnant may, un- 
der some circumstances, form a relatively cold disk 
that evolves without a substantial mass infusion at 
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its outer edge. Otherwise, there might also be a fos- 
sil accretion disk around the central black hole which 
is the remnant of a phase of higher activity in the 
past (Falcke & Heinrich 1994). In both cases the disk 
would inevitably interact with the infalling wind and 
thus change the matter inflow and the structure of the 
accretion disk. The resulting observational signature 
will then depend strongly on the ratio of wind to disk 
accretion rates, the viscosity in the disk, the specific 
angular momentum of the infalling matter and the 
time scales involved. 

The broad-band radiative emission from Sgr A* 
may be produced either in the quasi-spherical accre- 
tion portion of the inflow (Melia 1992 & 1994) with 
M w ~ i o 21 ~ 22 g s _1 , or by a combination of disk plus 
radio-jet (Falcke et al. 1993a&b; see Falcke 1996 for 
a review). As discussed above, quasi-spherical accre- 
tion seems to be unavoidable at large radii, but the 
low actual luminosity of Sgr A* seems to point to a 
much lower accretion rate in a starving disk and hence 
Sgr A* appears to be the ideal case for applying the 
concept of a fossil disk fed by quasi-spherical accre- 
tion. We note here that the disk may also be advec- 
tive, in which case the luminosity may be decreased 
substantially due to the large fraction of dissipated 
energy carried inwards through the event horizon by 
the disk plasma (e.g., Narayan ct al. 1995). In fu- 
ture developments of our model, we shall also consider 
such disk configurations. 

The outline of the present paper is as follows: in 
Sec. 2 we discuss the wind energy deposition onto the 
disk and we will derive the equations for the evolu- 
tion of an a-disk with arbitrary wind infall in Sec. 3; 
in Sec. 4 we will present numerical solutions for var- 
ious parameters and discuss their implication for the 
Galactic Center in Section 5. The paper is summa- 
rized in Sec. 6, where we also describe further appli- 
cations of our model. 

2. Wind Energy Deposition Onto The Disk 

When deriving the modified disk equations below, 
we will need to understand, at least crudely, the ef- 
fects of the impacting wind on the disk, and where the 
energy is deposited. Since the inflow is usually highly 
supersonic near the surface of the disk (e.g., Coker et 
al. 1996, Paper II), it is reasonable to assume that the 
wind termination proceeds via a shock. To estimate 
the relevant time scales, we can make the rough ap- 
proximation that an axisymmetric wind approaches 



the flat disk (of radius i?d = 10 16 cm-i?i6) in free fall, 
with an accretion rate M w = 10~ 4 M Q • M_ 4 . The 
black hole mass is taken to be M, = 10 6 M Q • Mq. As 
we shall see, these parameters are typical for weakly 
accreting nuclei, such as Sgr A* at the Galactic Cen- 
ter. For a strong shock, the downstream flow has a 
velocity 
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cm 



sec 



and a density 



3 x 10 5 cm" 



M_ 



= , (2) 



irR 2 v sh m p y/M 6 Rl 6 
while the temperature of the shocked gas is given by 

T sh = ^=6xlO'K^. (3) 

KB -Kl6 

Behind the shock, the gas will have to settle onto 
the disk, which in our case is much cooler with a tem- 
perature Td of several 100 to 1000 K. In this situ- 
ation, the dominant radiative cooling mechanism of 
the hot shocked gas with energy density dE/dV = 
3/2fc B T s 

h^sh would be thermal Brcmsstrahlung (ebrems 
1.7 x 10~ 27 [r/K] 1 / 2 [n/cm" 3 ] 2 crgsec- 1 cm- 3 ), cor- 
responding to a cooling time scale of ibrcms = 3 x 
10 9 sec M§Riq/ M_4. By comparison, the subsonic 
flow time scale is ta OV! = Hi/v s h, where Hi is the scale 
height of the shock above the disk. Thus, because 



= 13 



M, 



3/2 



(4) 



and Hi -C R, the gas cannot radiate appreciably be- 
fore reaching the disk. If we express R in units of 
the gravitational radius, this ratio is inversely propor- 
tional to the Eddington luminosity of the black hole, 
and can become smaller than unity only for super- 
Eddington wind accretion or very large radii. In the 
post-shock region above the surface of the disk, the 
state of the gas is then almost always dictated solely 
by the influence of Coulomb scatterings and conduc- 
tive heat transport. 

Let the temperature, density, and height of this 
interface region between shock and disk be 7], m, and 
H u respectively, where T s h > 7] > X"d and presumably 
Hi R (see below). The density is given roughly 
by the balance of ram pressure and thermal (surface) 
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pressure: 



Hi 



2.8 x 10 9 cm" 



M 6 1/2 M_ 4 / Ti 



p5/2 VlOOOK 
n 16 



(5) 



The conductive heat flux in the interface region is q ~ 
2.5 x 1CT 6 erg/sec/cm 2 (T;/K) 7 / 2 / (Hi/cm) (Frank, 
King & Raine 1985, hereafter FKR, Eq. 3.42). It 
is not difficult to show that the ratio of the thermal 
conductivity time scale i con = H ■ (dE/dV)/q to the 
flow time i; = H/v^ is 



M 6 5/2 



.3/2 



(6) 



which in most cases is much smaller than unity if 
Hi <C R, remains valid even for Eddingtion accretion 
and large radii, and is independent of the disk tem- 
perature. 

An upper limit to the value of the scale height 
Hi is given by the Coulomb energy exchange length, 
based on the scattering time scale i cou i = 4.6 x 
lO-^t^n^Tj^Kse^cm- 8 (FKR, Eq. 3.32) on 
which incoming particles share their energy with the 
background plasma. Since this time scale is extremely 
short, the wind/disk interface region is extremely thin 
(even compared to the disk scale height): 



f^2xl0- 



M, 
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(7) 



Thus, because of the high efficiency of thermal con- 
duction and Coulomb scatterings, and the relatively 
large ram pressure, in almost all cases of interest most 
of the energy in the flow is deposited directly in a very 
thin layer onto the surface of the disk. Since the sur- 
face of the disk (in contrast to its interior) may be only 
weakly ionized, the incoming ions might actually pen- 
etrate somewhat below the disk surface until Coulomb 
interactions start to work and as Hi is so small, most 
of the energy will be deposited somewhere in or even 
beneath the surface layer. Given these conditions, we 
shall treat the disk as a large mass and heat sink for 
the inflowing wind. In particular, we shall assume 
for simplicity that the deposited energy is radiated as 
a (local) black body, since the spectral modifications 
due to bremsstrahlung emission above this region are 
expected to be small (see Eq. 0). 



3.1. Mass and angular momentum conserva- 
tion 

The basic equations describing the evolution of 
an accretion disk are the mass conservation equation 
and the angular momentum equation (von Weizsacker 
1948; Lynden-Bell & Pringle 1974; see also Equa- 
tions 5.3 and 5.4 in FKR). In the case of an accre- 
tion disk intercepting an external wind, these have to 
be modifed to include additional source terms from 
the external gas flow. If the wind produces a mass 
infall E w (r,t) per unit area and time with a velocity 
u = (u r ,M0,u z ) in cylindrical coordinates, we obtain 
the mass equation 



and the angular momentum equation 



(8) 
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dr 
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The other variables used here are the disk surface den- 
sity E, the disk radius r, time t, (local) wind angu- 
lar momentum flux 1= rE w u^, disk internal viscous 
torque G=27rrj/Er 2 a/, partial radial derivative ' = 
d/dr, internal disk viscosity v, and the disk angular 
velocity uj=y // G gTav M/r 3 . Combining these equations 
yields 



dt 
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r dr 
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dr V 
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1/2 



o(l-0" 



(10) 

which is the well known partial differential equation 
for the time evolution of the surface density of the 
disk, but now including the effects of wind infall. The 
second equation, obtained by combining Equations |^ 
& 0, is the radial velocity of the accreting gas 
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rE 



dr 



(ii) 

To permit a straightforward identification of the 
physical meaning of the new terms, and to provide 
a means of comparing the various effects in a simple 
manner, we define an external, rotational 'viscosity' 
v w that describes the braking of the disk by the wind: 



3. Modified Disk Equations 



(i-C)V^, 



(12) 
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and a parameter £(r, t) that describes the angular 
velocity of the wind relative to its Keplerian value 
Vcf, = raj, such that 

U0 = ?^ . (13) 

For the type of accreted wind discussed in this paper, 
we have — 1 < £ < 1. The term (1— £) 2 in Equation |l2| 
ensures that the viscosity is always positive even if the 
angular velocity of the external wind is larger than 
the Keplerian value. In the case of a wind flowing 
out from the center intercepting the disk with £ > 1, 
the disk would be spun up and the radial disk inflow 
would be diminished by the external wind. This is 
reflected by a negative forefactor (1 — £) inside the 
brackets of Equations [h] & Ell The only disadvantage 
of this definition of u w is the apparent singularity at 
£ = 1, which, however, is not real as then also v w = 0. 

In this definition we only consider the viscosity in- 
duced by the azimuthal velocity of the wind. The 
radial contributions of the momentum equation are 
neglected here but are considered in the energy equa- 
tion as discussed below. 

3.2. Energy equation 

Any kind of viscosity, which is needed for the an- 
gular momentum and mass transport in the disk, is 
intimately linked to friction and dissipative processes. 
In accretion disk theory, it is usually assumed that the 
dissipated energy is radiated locally, though advection 
can under certain circumstances carry a significant 
fraction of the flux. The rate of energy dissipation 
is a function of the disk surface density, the viscosity 
and the shear rate. For the internal energy dissipation 
caused by turbulent viscosity and differential rotation 
in the disk, the rate of dissipated energy per unit area 
and time is (see FKR, Equation 4.25) 

D d = ^vZ(ru;') 2 . (14) 

This is defined as dissipated energy per surface area 
and as we have 2 sides of the disk we get a factor 
of 1/2. Braking of the disk by the wind will cause 
a change in angular momentum, lead to friction be- 
tween wind and disk and finally to energy dissipa- 
tion. For now we will neglect the vertical structure of 
the disk and wind and treat the disk-wind interaction 
zone as an infinitesimally thin sheet in which the en- 
ergy is dissipated (see Sec. 2 above). The dissipation 
rate, which is the change of energy AE per time At 



and surface area AA due to the friction between the 
wind and disk, can then be written as 

1 AE 
w " 2AtAA 

1 2ttt{AG w /t) 

2 At AA 

= g ( v 4> - u 4>) 2 ^™ i ( 15 ) 

using the torque 

AG w ^r± w (v 4> -u 4> )AA (16) 

exerted by the infalling wind on the accretion disk on 
a time scale 

At = 2irr/(v 4> -U4 > ) , (17) 

given by the azimuthal velocity of the wind relative 
to the disk. Thus we obtain the final form for the 
dissipation rate and can express this in terms of the 
rotational viscosity v w defined above: 

D w = (l/2)rV(l-0 2 (Sw/S)S 

= (l/2)i/ w Ew 2 . (18) 

Comparison of Equations Il8| and [l4| shows that our 
definition of v w (Equation jl2|) is indeed justified. 

In addition to the dissipative heating by rotational 
and turbulent friction we have to take the ram pres- 
sure of the wind into account, which results in a heat- 
ing rate of a form similar to Equation [l8| We shall 
consider both the u r and u z components of the wind 
ramming into the disk. It is assumed that the radial 
velocity of the disk is much smaller than the wind 
velocity v r <C yzif + uf and that the radial pressure 
gradient induced by the wind is negligible for the disk 
structure. If all the kinetic energy associated with u r 
and u z is converted into heat and thus contributes to 
the total dissipation rate, we have an additional term 

D lam = ^± w (t> 2 + vl) . (19) 

The total dissipation rate per surface area is then 
given by 

D = D d + D w + D lam . (20) 

One should note that in any realistic case, the wind 
properties will be asymmetric with respect to the two 
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sides of the disk: one would have to split Equations [l8| 
and |l9| accordingly and introduce separate parameters 
£i|2, S wl | 2 and v w \\2 an d obtain asymmetric dissipa- 
tion rates. We do not want to introduce this compli- 
cation into our model at the current stage. 

3.3. The vertical structure 

To solve the disk equations and obtain the evo- 
lution of £ we need to know the radial dependence 
of the turbulent viscosity v. This is usually done by 
introducing the Shakura & Sunyaev (1973) a param- 
eter, defined such that 



v = ac s z 



(21) 



where c s is the local sound speed, zo is the scale height 
of the disk and a is the dimensionless viscosity param- 
eter, usually between and 1. This has the benefit 
that the viscosity can be expressed in terms of local 
properties of the disk but has the disadvantage that 
one needs to solve its vertical structure. 

The disk scale height can be obtained by demand- 
ing pressure equilibrium between internal gas pres- 
sure, gravitational attraction in the z-direction and 
external ram pressure from the wind. Substituting 
the gradient d/dz with the scale height zq, we obtain 
for a gas-pressure dominated disk with sound speed 



\ 




1 - 



1 2j w U z 

2 YiOj c s 



(22) 



This means that the disk is substantially squeezed if 
u z S] w /5] 3> 2c s w; if u z is of the order of the Keplerian 
velocity (i.e., u z ~ roj) then we have the condition 



2£/£ w < (c s r)" 



(23) 



for the time scale involved in order to compress the 
disk. The left hand side is the time scale of the wind 
interception and the right hand side the global thermal 
time scale of the disk. 

To complete our set of basic equations, we need to 
specify the central temperature in the disk to get the 
sound speed. We assume that the energy dissipated 
by viscous processes in the disk diffuses outwards from 
the central plane by radiative energy transport in an 
optically thick medium (with optical depth r 3> 1) in 



local thermal equilibrium. This will give us an energy 
flux of dissipated energy in the disk at its surface of 



F d (z Q ) = D d = 



4a AT(z ) 4 



(24) 



3 At 1 

the A here signifying changes in the corresponding 
quantity between the disk's midplane and its surface. 
The total energy flux will then be the sum of the en- 
ergy dissipated in the disk and the energy dissipated 
by the impact of the wind on the surface of the disk. If 
this is radiated immediately by black body radiation 
with an effective temperature T e g, we have 

F tot (z ) = aT cff (r,t) 4 = D(r,t), (25) 

and we can set T(z ) = T cS . With AT = T(z ) - 
T(0) and Ar — r, we can solve Equation ||| for T(0), 
getting 



T(0) 4 



3r 

T 



1 Da 



D ra 



(26) 



This equation reduces to the usual equation for the 
central temperature if D w ~ 0, D xam ~ and t>1, 
but it also allows a situation where the dissipation 
is dominated completely by the impact of the wind. 
This will increase the central temperature and make 
the temperature gradient shallower and in the ex- 
treme case could lead to a quasi-isothermal disk; in 
the cases discussed in this paper, this does not hap- 
pen. 

With the central temperature at hand we can now 
evaluate the pressure 



P 



and the sound speed 



fim p 



(27) 



(28) 



if we solve the second order partial differential equa- 
tion for dTi/dt (Equation [l0[). We only need to specify 
the opacity. Assuming that we are in the outer and 
cooler region of the disk, we may consider simply the 
contribution from free-free absorption which gives us 



6.6 ■ lO^cmVg 




-7/2 

T^7) • (29) 
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This opacity will become progressively less accurate 
as we move to very cold or very hot disks, but we 
ignore this trend here, since we are primarily inter- 
ested in a basic understanding of the disk evolution 
and the main effects. We hope to include a more elab- 
orate opacity function for more realistic calculations 
in a future version. 

The complete set of equations can now be solved 
numerically. For each time step we calculate the ra- 
dial gradients in £ on the right-hand side of Equa- 
tion to get dT,/dt and an estimate for £ at the 
next time step using a simple Runge-Kutta algorithm 
with adaptive step size. The vertical structure has to 
be calculated at every step and because the height is 
no longer a simple function of c s and lj, we have to 
solve this set of equations implicitly. Fortunately, c s 
changes only slowly so one can use the sound speed 
from the previous step as a first guess and achieve 
convergence very quickly. 

3.4. Spectrum 

To get a rough spectrum from the accretion disk, 
we have integrated black body spectra with the local 
effective temperature T c ff (r) over the entire disk sur- 
face. Following our discussion in Sec. 2, we consider 
T e g to represent the total dissipation rate, as indi- 
cated in Equations & Here, we have neglected 
all vertical transport processes that might modify the 
disk spectrum. Mainly because of the difficulties in 
comparing disk spectra to observations and the uncer- 
tainties involved in the vertical disk structure itself, 
the theory for radiation transport in disk atmospheres 
is not well developed, especially in the case of AGN 
disks. While Sun & Malkan (1989) find that the pure 
black body integration can explain AGN spectra suffi- 
ciently well, Laor & Netzer (1989) have used a simpli- 
fied treatment to account for the effects of the vertical 
structure. The changes to the calculated spectrum 
are, however, so small, especially compared to the 
scarce data for the Galactic Center, that inclusion of 
these effects would neither be opportune nor testable. 
For simplicity we also calculate our spectra for only 
one inclination angle i = 60° (i.e. cosi = 0.5). The 
spectra we show should therefore be considered only 
as a first order approximation and one should bear 
in mind that orientation effects, relativistic effects, 
limb-darkening, and propagation effects will lead to a 
somewhat different appearace. 



4. Numerical Solutions 

As an example of the evolution of a disk with wind 
infall, we choose parameters that might roughly be 
appropriate for the Galactic Center, i.e. a black hole 
mass of M. = 10 6 Mq, an accretion rate for the fossil 
disk of Mq = 10~ 8 MQ/yr, a wind infall rate of M w = 
lO _4 M0/yr, and an a of 10 -4 (e.g. Falcke & Heinrich 
1994). For simplicity in this first series of tests, we 
will correlate the vertical wind velocity directly to its 
radial velocity, such that 



u z = y/l — £ 2 «0 and u T = 0. (30) 

The boundary conditions for the innermost ring of 
the disk (n n = QR g = 6G gra v-/Vf./c 2 ) are such that 
we impose dH/dt = and £ = to account for 
the presence of the black hole. The outer boundary 
(at r = r ou t) can be free (using only the one-sided 
derivative) or be fixed to a certain M(r out ) by virtue 
of Equation [□] and M(r) = — 27rru r £, thus setting 
the square bracket in Equation [lO] to M (r out ) /6tt at 
r = r out . The radial velocity is also limited to not 
exceed the free-fall velocity. Even though we are here 
discussing the situation pertaining to a black hole, we 
ignore any relativistic effects, which are not critical as 
long as the wind affects mainly the outer portions of 
the disk. During the calculations we check that the 
total mass in the system is indeed conserved. 

4.1. Validation of the code 

To validate our code and test its reliability, we first 
calculate the disk evolution in two scenarios without 
wind infall to compare the surface density profiles 
with those of the stationary solution. 

Empty Disk We start with an almost empty, very 
low £ disk and feed it with a very high accretion rate 
close to the Eddington value of M = 10~ 2 MQ/yr 
at r out until an equilibrium state is reached. Figure [j] 
shows how the accreted matter quickly fills the almost 
empty disk and then slowly approaches the analytical 
Shakura & Sunyaev solution for the steady state (up- 
per panel in Fig. [j]). Also, M approaches a constant 
value for all r. 

Gaussian Ring In a second test, a textbook ex- 
ample, we start with a ring with a Gaussian surface 
density distribution in r and no accretion from the 
outer boundary. Figure shows that the ring evolves 
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towards the stationary solution with constant M. In 
the outer regions, matter is actually flowing out, car- 
rying away the angular momentum; only there one 
can see a deviation from the stationary solution. In 
the later evolution (not shown here), the disk runs 
out of fuel with a decreasing surface density and size 
but it retains its basic shape. 

The agreement between the analytical and numer- 
ical solutions demonstrates that our code produces 
reliable results. 

4.2. Disk evolution with a large scale wind 

In the first scenario where we allow a wind to fall 
onto the disk we assume that the disk intercepts the 
inflow over a fairly large region, extending out to a 
radius r — 90000i? g ~ 10 16 cm. In the absence of 
detailed hydrodynamical simulations, we do not yet 
know what the radial distribution of the wind is and 
we here consider the case of a constant mass infall 
per unit area, which means that the impact of the 
wind is felt most strongly at the largest radii, but 
is still present in the inner regions. More realistic 
distributions will be discussed in the second paper, 
where we will use the results of numerical 3D wind 
calculations. 

Maximum Angular Momentum Figure || shows 
the evolution of a fossil disk intercepting a wind with 
constant E w over a large collecting described 
in the previous paragraph. Hence, the largest frac- 
tion of the wind is intercepted at large radii, which 
leads to a significant increase of the surface density 
there. Because the viscous time scale (starting with 
the assumed a = 10~ 4 ) is too long for a radial trans- 
port of the matter — the viscosity increases only by 
a factor of three during the simulation — the modest 
increase in £ at smaller radii is due mostly to the lo- 
cal wind infall. After many time steps, the pile-up in 
the outer disk starts to become less pronounced due 
to radial accretion and the entire disk begins to con- 
verge towards the stationary solution. The time scale 
for an increase in surface density and luminosity is 
larger than about 10 6 yrs. We note that the dynam- 
ical time scale (~ 20 years) for the disk calculated 
from its rotation velocity is 10 5 times shorter than 
the time scale for the disk structure to change in re- 
sponse to the wind. Hence, the most notable change 
is expected to be the spectrum. The wind infall leads 
to a local increase in the disk accretion rate and the 



production of an infrared (IR) bump — the spectrum 
becomes almost flat. 

Hurricane If we simulate a "hurricane" (Fig. [|), 
where all the mass is deposited within an annulus 
(bound by the radii 30000-90000 R s ), the mass piles 
up in the outer region with a gradual transfer toward 
smaller radii over a 10 8 -yr time scale. Again, this re- 
flects the large viscous time scale in the disk assuming 
that a = 10~ 4 remains constant over this period, and 
is similar to the 'empty disk' scenario discussed above. 

Zero Angular Momentum The time scale for the 
transfer of mass inwards toward smaller radii short- 
ens dramatically when the infalling wind carries sub- 
Keplerian angular momentum. The importance of the 
variations in angular momentum and velocity distri- 
butions in the wind is highlighted in Figure |5|, where 
we have repeated the previous calculation (shown in 
Fig. 3) for a free-falling wind with zero angular mo- 
mentum (£ = 0). The IR bump jumps almost im- 
mediately to the value that corresponds to the high 
accretion rate in the wind. The reason for this is that 
the non-zero vertical impact velocity of the wind and 
the strong friction between wind and disk lead to a 
rapid dissipation of energy and rapid radial accretion. 
The deposition of mass with zero angular momentum 
in the disk also leads to a compression of the disk in 
the radial direction (the outer rings move inwards as 
they become sub-Keplerian) and results in an increase 
in surface density beyond the value induced by the in- 
creased mass deposition through the wind alone. Also 
the viscosity increases much more strongly. 

4.3. Disk evolution with a small scale wind 

In some systems, the wind flowing into the inner 
disk region may carry very little specific angular mo- 
mentum. In this case, because the circularization ra- 
dius can be significantly smaller than the fossil disk's 
outer radius, the impact of the infalling gas on some 
portions of the disk can be much larger because the 
wind is now restricted to a smaller disk surface area 
and hence accounts for a significantly larger local flow 
rate. We have therefore repeated the above calcu- 
lations, but now with a reduced interception region 
designated by a radius 100 times smaller than used in 
the previous section. This leads to a much stronger 
effect on the disk because of the higher £ w and the 
correspondingly shorter time scales. 
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Maximum Angular Momentum For a wind with 
the local Keplerian velocity (Fig. ||), the disk reacts 
almost immediately at all radii because the wind in- 
fall dominates the radial accretion rate everywhere. 
Still, the infall is more prominent in the outer region 
in the beginning, yielding an almost radially constant 
surface density. As soon as the viscous mass transport 
has caught up, the surface density profile approaches 
the stationary solution. The mass transport at the 
outer wind radius is bi-directio nal , like in the gaus- 
sian density distribution (Sec. 4.1), so that beyond 



the wind accretion radius matter is flowing outwards 
and not inwards. The spectrum is immediately dom- 
inated by emission from the outer region because the 
kinetic energy of the wind (in u z ) is dissipated im- 
mediately, but since the interaction region of the disk 
is so small, we do not see two bumps. After viscous 
processes have transported the higher accretion rate 
inwards, the spectrum simply broadens somewhat to- 
wards higher frequencies. 

Zero Angular Momentum The situation becomes 
even more extreme if the wind has no angular momen- 
tum. Within ~ 1000 years (which is still much longer 
than the dynamical time scale) the incorporated mass 
in the outer region becomes equal to the stored mass 
in the disk. Because of the low angular momentum, 
the mass falls into a lower orbit within this short pe- 
riod of time, leaving a gap behind until the fossilized 
angular momentum in the disk stops a further free 
fall. The accretion flow is now purely inwards and af- 
ter it is stabilized continues to grow on a longer time 
scale (10 5 yrs). The spectrum is immediately dom- 
inated by a single Planck-spectrum from the outer 
region. 

5. Preliminary Application To Sgr A* 
5.1. The basic scenario 

As our first (and necessarily simplified) applica- 
tion of the results in the previous section, we now 
consider the situation in the Galactic Center. There 
exists ample evidence for the presence of a strong 
wind in the region around Sgr A*, with a charac- 
teristic velocity v w ~ 500 — 700 km s -1 and mass 
loss rate M w ~ 10~ 3 M Q yr" 1 (see Melia 1994 and 
the references cited therein, particularly Hall, Klcin- 
mann & Scoville 1982, Geballe et al. 1991, Krabbe 
et al. 1991 and Gatley et al. 1986, among oth- 
ers). In addition, the latest stellar kinematic studies 



point to a likely mass of M « 1.8 (±0.5) x 10 6 M Q 
for the central object (Haller et al. 1995), and so 
it is expected that Sgr A* is accreting at a rate of 
about 10(- 5 ) to (- 4 )M Q /yr (Ruffert & Melia 1994) 
via Bondi-Hoyle accretion. 

However, the luminosity of Sgr A* is tightly con- 
strained to a value of a few 10 5 Lq (Falcke et al. 
1993a; Zylka et al. 1995). Thus, if the accretion flow 
eventually turns into a disk at smaller radii, and as- 
suming a conversion efficiency of about 10% rest mass 
energy into radiation, we infer a disk accretion rate 
of a few 10 _7 MQ/yr. If the winds surrounding the 
central black hole are relatively uniform, the fluctua- 
tions in the accreted specific angular momentum are 
sufficiently small that only a small (less than about 
5 — 10 Schwarzschild radii in size) disk is then ex- 
pected to form (Melia 1994; Ruffert & Melia 1994). 
Such a small disk could be consistent with the low 
luminosity observed from this region, but it is likely 
that gradients in the ambient medium produce a big- 
ger specific angular momentum in the accreted mate- 
rial, and hence a larger circularization radius. This 
could be due in part to the fact that the Galactic 
Center wind is comprised of several stellar compo- 
nents, which would make the flow around Sgr A* non- 
uniform. Departures from uniformity increase the ac- 
creted specific angular momentum since the cancella- 
tions in the post bow-shock flow are not complete — 
one side tends to dominate over the other, producing 
a net accreted angular momentum. The correspond- 
ingly larger circularization radius produces a bigger 
and presumably brighter disk. Another way of say- 
ing this is that the infalling gas from a non-uniform 
medium retains a larger Keplerian energy that must 
be dissipated as it drifts inwards toward the black 
hole. 

It could be that the disk, if large, is highly ad- 
vective (Narayan et al. 1995). The role of advection 
has usually been ignored in a-disk theory, but sev- 
eral studies (e.g., by Abramowitz et al. 1988, 1995, 
and more recently by Narayan & Yi 1994) have shown 
that under some circumstances, the internal disk en- 
ergy can constitute a substantial radial flux when it 
is transported by the inward drift. As much as 90% 
or more of the locally dissipated gravitational energy 
flux can be carried away in this manner, thereby de- 
creasing the radiated luminosity for a given mass ac- 
cretion rate by as much as an order of magnitude, and 
possibly more. However, without detailed radiative- 
hydrodynamical simulations, this scenario is at best 
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still only a possibility. 

It has also been suggested that Sgr A* might 
be surrounded by a large, slowly accreting, fossil 
disk, perhaps formed from the remnants of a tidally- 
disrupted star that ventured too close to the black 
hole. In this case one expects that the large-scale 
Bondi-Hoylc inflow must then be captured and in- 
corporated into the disk. The latter might remain 
faint if most of the quasi-spherical infall circularizes 
at large radii, for the incorporated gas would require 
a fairly long viscous time scale in order to reach the 
central black hole and produce the correspondingly 
high luminosity output (Falcke & Heinrich 1994). A 
fossil disk could also act as a resevoir for the infalling 
plasma, thereby delaying the dissipation of its energy. 

5.2. Model vs. observation 

Many of the calculations reported in this paper cor- 
respond to a situation similar to that of Sgr A*. In 
order for a particular model to be relevant to this 
source, it should be able to account for the high 
wind-accretion rate and yet a low luminosity of a few 
10 5 L Q . Recent NIR observations of Sgr A* appear to 
show 5 different sources within the Sgr A* error box 
(Eckart et al. 1995), though the source separation 
in this image is smaller than the telescope's diffrac- 
tion limit. If this separation of the infrared emission 
from this region is real, the limit on Sgr A*'s lumi- 
nosity could be as low as L u < 3 • 10 21 erg/sec/Hz 
(for an interstellar absorption of Ak — 3.7) in the 
K-band (~ 10 14 Hz). The upper limit for Sgr A* at 
12.4/im (10 13 - 4 Hz) is 0.1 Jy (0.9 • 10 22 erg/sec/Hz) 
(Gezari et al. 1996). Recently Stolovy, Hayward, & 
Herter (1996) claimed a detection of Sgr A* at 8.7/xm 
(100 mjy dereddened) which may be interesting with 
respect to our work but needs further confirmation. 
A starving accretion disk with M < 10 -8 Mq/yt, 
as used in our simulations, could be consistent with 
those numbers when the following additional factors 
are taken into account. Requiring that the disk is 
seen edge-on (cos 8 <; 0.1), would reduce the flux dis- 
played in the figures by another factor 5. Moreover, 
there are indications that Sgr A* is intrinsically ab- 
sorbed (Predehl & Triimper 1994) — in addition to 
the interstellar attenuation — which would reduce the 
observed luminosity further. 

A second important factor to take into account 
when comparing particular models to the observa- 
tions is the time scale required for a change in the 
radial accretion rate. If the main components in the 



Bondi-Hoyle accretion are the stellar winds of cen- 
tral He I stars, the accreting medium should have 
been there a certain fraction of the stars' lifetime 
(~ 3 - 7 x 10 6 years, Krabbe et al. 1995), though 
probably not longer than this. We would therefore 
expect that a typical time scale for the accretion pro- 
cess should be longer than of order a million years if 
a fossil disk is involved. 

The luminosity and time scale constraints would 
therefore seem to rule out situations in which the wind 
is captured by the disk at small radii, because they 
produce spectra that are too luminous in the NIR 
and IR very quickly after the onset of Bondi-Hoyle 
accretion (Figs. |^ & 0). In other words, the post-bow 
shock flow cannot be carrying a small specific angu- 
lar momentum (which would lead to circularization 
at small radii), but rather should have sufficient az- 
imuthal velocity to allow it to merge with the fossil 
disk near its outer edge. However, models with a low 
angular momentum wind that is captured far out in 
the disk fare no better (as far as Sgr A* is concerned; 
see Fig. ||) because they violate the IR limits due 
to their instant transformation of kinetic energy into 
heat and thence into IR radiation. 

It appears from our analysis of the calculations 
that the only class of models that might be consistent 
with the observations of Sgr A* are those in which the 
angular momentum of the wind is large enough that 
the infall circularizes at a radius > 10 16 cm and set- 
tles gently onto the disk (Figs. H°& ||). ln addition, 
the a-parameter must be small (a <C 1) in order to 
prolong the transfer of the deposited wind material 
into the center. The 'hurricane' model with a <J 10 -4 
seems to be a viable scenario, though this model does 
not even require the presence of a fossil accretion disk 
to stop the infall because the high angular momentum 
wind will circularize and stop at larger radii anyway. 
If it exists, a starving, fossil disk may even lie closer 
to the black hole and not be coupled directly with the 
large scale disk formed by the wind infall. 

6. Discussion 

In this paper we have rederived the basic equa- 
tions for an a-accretion disk with the consideration of 
wind infall. Besides a source term for the mass, the 
wind infall also introduces a dynamical viscosity v w 
to the basic equations when its angular momentum 
is sub-Keplerian. We have calculated the evolution 
of an accretion disk with wind infall for various sets 
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of parameters, leading to some interesting effects not 
present in conventional accretion disk theory. The 
situation we considered in detail is that of a fossil 
accretion disk, which intercepts a wind with a much 
higher accretion rate than the former. As long as 
the viscosity parameter in such a disk is very small 
(a <C 1), the time during which the surface density of 
the disk changes is much longer than the dynamical 
time scale because of its large mass storage capability. 
In this case, the disk is not blown away and it can in- 
corporate winds with very high mass accretion rates. 
A large, non-coaxial angular momentum in the wind 
may, however, lead to a warping of the disk which we 
ignore here. 

The overall evolution of the disk may be described 
in terms of two basic classical solutions: (a) that of 
an almost empty disk that is fed from the outside and 
(b) an initial ring of matter that evolves with time 
and where most of the matter is flowing in, while a 
small fraction flows outwards carrying away the an- 
gular momentum. If the intercepted wind has a low 
specific angular momentum, the disk will shrink in 
the radial direction on a short time scale determined 
by S/S w and de-couple from its outer region, form- 
ing a gap. The matter piles up until viscous processes 
in the disk become important and the inner regions 
are filled up like in the empty disk evolution. For a 
high specific angular momentum wind, the disk will 
also have a radial outflow in its outer portion like in 
the ring evolution. The interaction between wind and 
disk leads to an additional heating of the disk on its 
surface and a moderate increase in the viscosity and 
therefore a faster evolution; the effect is stronger for 
a sub-Keplerian wind. 

The most important changes with respect to the 
conventional solutions concern the emitted spectrum: 
the emission at lower-frequencies is enhanced. This 
may produce either a flattening of the spectrum or 
even a strong low-frequency bump (in AGN the bump 
would be produced in the IR). This bump can be as 
luminous as the primary bump produced in the inner 
region of the disk or even brighter. The most drastic 
effect is due to the angular momentum deposited by 
the wind. When the wind has a Kcplcrian velocity 
and settles gently onto the disk (it z = 0), the forma- 
tion of the low-frequency bump may take a long time 
S/E w , while a wind with a low specific angular mo- 
mentum and a high u z will release its kinetic energy 
very quickly and the low-frequency bump will appear 
almost instantaneously. In both cases, however, it 



will take a viscous time R 2 jv for the increased accre- 
tion rate in the outer part to reach the inner region 
and lead to an increase of the (usually dominating) 
high-frequency emission. 

The process described here has many potential ap- 
plications. The production of an IR bump by an in- 
falling wind, perhaps from the nuclear star cluster or 
a circumnuclear starburst, would be an interesting 
although as yet undeveloped — signature to inves- 
tigate. A scenario such as this can also be invoked 
to account for some observational characteristics in 
winds onto protostellar disks. In this paper, we have 
briefly discussed the application of our model to the 
Galactic Center source Sgr A*. 

The presence of strong stellar winds and a high 
central mass concentration in the Galactic Center ap- 
pear to be inconsistent with the low Infrared luminos- 
ity of the putative supermassive black hole, Sgr A*. 
We have considered the possibility that this object 
may be surrounded by a fossil disk accreting at a rate 
of 10 _8 MQ/yr, while it is intercepting the incoming 
wind with M w <J 10 _4 MQ/yr. Current observational 
limits for Sgr A* are very restrictive and we have 
found that they are not consistent with most of the 
wind specific angular momentum configurations we 
have considered here. As the estimates for the bolo- 
metric luminosity of Sgr A* may change, the direct 
observations in the IR and NIR are the most critical 
limits to be observed in this respect and we find that 
any fossil disk would become too luminous too quickly 
in those bands. We can only reconcile the data with 
the models if the wind circularizes on a scale of ;> 10 16 
cm and if the a-parameter within the disk is <; 10~ 4 . 
Without detailed numerical simulations, these num- 
bers are at best only estimates. Even so, it seems 
safe to conclude that if a fossil disk is present around 
Sgr A*, then it must have a very low a and the wind 
must carry a large specific angular momentum in or- 
der to describe Sgr A* self-consistently. This rules out 
a large region of the phase space available for possible 
models of Sgr A*. 

In a future extension of this work, we plan to couple 
our calculations to 3D hydrodynamical simulations 
for the large scale Bondi-Hoyle accretion (Paper II) 
in order to get a better estimate for the radial distri- 
bution of the wind infall. As we have seen here, the 
specific angular momentum of the wind appears to be 
the most critical parameter. Further improvements of 
our model will include a more accurate description of 
the opacity and the possible warping of the disk by a 
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very strong, asymmetric wind. 
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Fig. 1. — Surface density evolution of an almost 
empty disk that is fed from the outer boundary with 
lO~ 2 M /yr. The analytical stationary solution is 
given by the dashcd-dottcd line. The lower panel 
(the zero £ distribution is outside the plot range here) 
shows the evolution on an extended time scale com- 
pared to the upper panel. 

Fig. 2. — Surface density and accretion rate evolu- 
tion of a Gaussian ring with no external feeding. For 
comparison, two appropriately scaled stationary so- 
lutions are plotted, the upper curve corresponding to 
M = 10~ 2 M Q /yr and the lower one to an accre- 
tion rate roughly 5 times lower, which is the actual 
accretion rate in the disk. 

Fig. 3. — Surface density and emitted spectrum 
evolution of a fossil accretion disk with strong wind 
deposition; the wind has Kcplerian azimuthal veloc- 
ity (M. = 10 6 M Q , M = l(T 8 M /yr, a = 1(T 4 , 
M w = lCT 4 M Q /yr, £ = 1, and r w = 90000i? g ). 

Fig. 4. — Surface density and emitted spectrum evo- 
lution of a fossil accretion disk with strong wind de- 
position only within the outer ring (a "hurricane"); 
the wind has Keplerian azimuthal velocity (M. = 
1O 6 M , M = 10- 8 M o /yr, a = 1CT 4 , M w = 
lCT 4 M Q /yr, £ = 1, and r w = 30000 - 90000i? g ). 

Fig. 5. — Surface density and emitted spectrum 
evolution of a fossil accretion disk with strong wind 
deposition; the wind has zero azimuthal velocity 
(M. = 10 6 M Q , M = 10- 8 M Q /yr, a = 10~ 4 , 
M w = 10- 4 M Q /yr, £ = 0, and r w = 90000i? g ). 

Fig. 6. — Surface density, emitted spectrum, and 
accretion rate evolution of a fossil accretion disk with 
strong wind deposition at small radii; the wind has a 
Kcplerian azimuthal velocity (M, = 1O 6 M , Mo = 
10- 8 M Q /yr, a = lO" 4 , M w = lO- 4 M /yr, £ = 1, 
and r w = 900i? g ). For comparison we have plotted 
a stationary solution on top of the S curves (dashed- 
dotted line, top figure). 

Fig. 7. — Surface density, emitted spectrum, and 
accretion rate evolution of a fossil accretion disk 
with strong wind deposition at small radii; the wind 
has zero angular momentum (M. = 1O 6 M , M = 
10- 8 M Q /yr, a = 10" 4 , M w = lO- 4 M /yr, £ = 0, 
and r w = 900i? g ). For comparison we have plotted 
a stationary solution on top of the S curves (dashed- 



dottcd line, top figure). The height in the 'gap' is an 
artifact, as there is no disk anymore in this region. 
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